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Abstract 



Parasupersymmetry of the one-dimensional time-dependent Schrodinger equ- 
Q\ ■ ation is established. It is intimately connected with a chain of the time-dependent 

Darboux transformations. As an example a parasupersymmetric model of non- 
Q_i! relativistic free particle with threefold degenerate discrete spectrum of an integral 

of motion is constructed. 

1. Supersymmetric quantum mechanics originally introduced by Witten |]J for in- 
vestigation of the supersymmetry breaking in quantum field theories attracts now a 
considerable attention from the different points of view (see a recent survey ||). Dif- 
ferent generalizations of the initial constructions are known. We can cite A-extended 
supersymmetric quantum mechanics Q, higher-derivative supersymmetry in quantum 
mechanics ||, and parasupersymmetric quantum mechanics 0-0. It is worth stress- 
ing that all these constructions deal with a stationary Schrodinger equation and con- 
sequently can be referred to the stationary supersymmetric quantum mechanics. The 
nonstationary supersymmetric quantum mechanics needs to be developed. First steps 
in this direction have been made in Refs. - H- 

An essential ingredient of the stationary supersymmetric quantum mechanics con- 



stitutes the well-known Darboux transformation [11] of the stationary Schrodinger 
equation in mathematics (see Ref. fTOfl). Recently the time-dependent Darboux trans- 
formation has been proposed |7j and on its base the supersymmetry of the nonstation- 
ary Schrodinger equation has been established J5]. This transformation seems to be 
very fruitful for construction of the nondispersive wave packets (coherent states) for 
anharmonic oscillator Hamiltonians with equidistant and quasiequidistant spectra ||. 

In this paper we establish the parasupersymmetry of the time-dependent Schrodin- 
ger equation which is intimately connected with a chain of the time dependent Darboux 
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transformations. With the help of these transformations we can construct a chain of the 
new exactly solvable nonstationary potentials starting from the initial potential even if 
it does not depend on time. We introduce the notion of completely reducible chain of the 
Darboux transformations, which is also valid for the stationary case. We associate with 
this chain a nonlinear algebra of the parasupersymmetric structure. As an example 
the application of two consecutive time-dependent Darboux transformations to the free 
particle Schrodinger equation (zero potential) is made. This leads us to two exactly 
solvable nonstationary potentials. For each potential the discrete spectrum basis of 
the Hilbert space defined over solutions of the nonstationary Schrodinger equation is 
obtained. Three nilpotent conserved supercharges together with the even symmetry 
operator formed a nonlinear algebra of the parasupersymmetric structure are explicitly 
constructed. 

2. In this section we summarize following Refs. 0-Q the main properties of 
the time-dependent Darboux transformation and introduce the supersymmetry of the 
nonstationary Schrodinger equation. 

Let us consider two time- dependent Schrodinger equations 

(id t - h ) V (x, t) = 0, h = -dl + V (x, t) , (1) 

(id t - h x ) <p (x, t) = 0, h 1 = -dl + V x (x, t) , (2) 

d 2 x = d x - d x , d x = d/dx, i El, x G R. 

The interval R = [a, b] for variable x can be either finite or infinite. Let T and T\ be lin- 
ear spaces of the solutions of Eqs. (JTJ) and (0) respectively defined over complex number 
field C. We assume that the solutions of Eq. ([J), called the initial Schrodinger equa- 
tion, are known. The problem the Darboux transformation solves consists in derivation 
of such a set of potentials V\ (x, t) that solutions ip (x, t) of Eq. (|j) could be obtained 
by means of acting of a first-order differential operator L called Darboux transforma- 
tion operator on the solutions ip(x,t): <p(x,t) = Lip(x,t). Since differential operators 
may have nontrivial kernels, we have the embedding T\ — {ip : ip — Lip,Vip G To} C T\. 
This signifies that Eq. (0) may have solutions which cannot be obtained by acting of 
operator L on any function ip G T . If necessary these solutions have to be obtained 
by other means. The condition that the operator L be first-order differential operator 
is very restrictive for the potential V\(x, t). We have to look for this potential together 
with the operator L. In order to achieve this goal (i.e. to find the unknown functional 
coefficients of the operator L and the function Vi(x,t)) we need to solve the operator 
equation H] 

L (id t - ho) = (id t - h) L 

in the space T . By virtue of Eq. (P the derivative of the function ip(x,i) G T 
with respect to t can be replaced by [—0% + V (x, t)} ip(x,t). Hence, if we will restrict 
ourselves to the first order differential operators we should not include the derivative d t 
in L and the general form of this operator should be as follows: L = a(x, t) + b(x, t)d x . 
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We refer the reader for details of calculations to Refs. |J and || and will give only 
the final result. The transformation operator L has the form 



L = L 1 (t)[-u x (x,t)/u(x,t) + d x ], (3) 

where the subscript x denotes the partial derivative with respect to the coordinate x 
and 



2 / d£lm(logw) a 



L\ (t) = exp 

The potential Vi(x,t) is defined by the potential difference 



(4) 



Ao t i(x,t) = h 1 -h = Vi_(x,t) - V (x,t) = -[logHx,^! 2 ]^. (5) 

Note that the operator L and the new potential Vi(x,t) are completely defined by 
a function u(x,t) called the transformation function. This function is a particular 
solution of the initial Schrodinger equation ([!]) subject to the condition 



(log 4) =«. ( 6 ) 

\ Us xxx 



U ' xxx 

called the reality condition of the new potential (the "*" means complex conjugation). 
Operator L + which is Laplace adjoint to L 

L + = -L 1 (t)[u:(x,t)/u*(x J t) + d x } (7) 

realizes the transformation in the inverse direction, i.e. the transformation from the 
solutions of Eq. (fj) to the ones of Eq. ([[]). As a result, the product L + L is a symmetry 
operator for Eq. (H) and LL + is the analogous one for Eq. (||). 

The action of the operator L (^) on the transformation function u gives zero. How- 
ever if this function is chosen in accordance with the Eq. (Q) we can readily verify that 
the function ^ 

V(x,t) — — t t r (8) 

is the solution of the new Schrodinger Eq. (Q). Note the following property of this 
function: L + v = 0. 

With the help of the transformation operators L and L + we build up the time- 
dependent nilpotent supercharge operators 

«-(2S). «-(o L o + ) w> 

which commute with the Schrodinger super-operator ild t — H where H = diag{/io, hi} 
is time-dependent super-Hamiltonian and I is the unit 2x2 matrix. In general, the 
super-Hamiltonian is not integral of motion for the quantum system guided by the 
matrix Schrodinger equation 

(ild t - H) * (x, t) = 0. (10) 
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The two-component function \1/ (x, t) belongs to the linear space Tqi defined over com- 
plex number field and for every if) G T Q spanned by the basis = ipe + , = Lipe_ 

rp rp 

with e + = (1,0) and e_ = (0, 1) , the superscript "T" stands for transposition. 

The operators (H) are integrals of motion for Eq. (|T0"|) . Using the symmetry op- 
erators L + L and LL + we can construct another integral of motion for this equation: 
S = diag {L + L, LL + }. The operators Q, Q + and S realize the well-known superalgebra 
[0 sZ(l /l) but instead of Hamiltonians we see in its construction another symmetry 
operators. In general, the operator S (as well as Q) depends on time and consequently 
we have the time-dependent superalgebra. 

3. The majority of cases of physical interest are such that it is possible to introduce 
the Hilbert space structure H°(R) in the space T with the appropriately defined scalar 
product. We suppose that the operator L + is the operator Hermitian conjugate to L 
with respect to this scalar product. The symmetry operator L + L for L being of the 
form (|3|) is of the second degree in d x . By means of Eq. (|TJ) the second derivative with 
respect to x for any if) G T can be replaced by the first derivative with respect to t. 
This means that the restriction of the operator L + L on the space T is an operator of 
the first degree in d x and d t and consequently the operator iL + L belongs to the real 
Lie algebra Go of symmetry operators for Eq. ([]]) which we assume to be known. Since 
the transformation function u is an eigenfunction of the operator iL+L corresponding 
to the zero eigenvalue we can denote L+L = ig(°> — a with g(°> G Go and ig^u = aou, 
ao G M. The operator g^ is usually skew- Hermitian and ig^ has a unique self-adjoint 
extension (see, for example, |12]]). In H°(R) this operator may have discrete spectra, 
continuous spectra or a combination of the two. The same also apply for the operator 
LL+ = igW — a 0) 9 ^n, G\ being a Lie algebra of symmetry operators for Eq. 
(□) . Operator gW is said to be a superpartner of g(°>. It is not difficult to calculate 
for the transformation function u subject to the reality condition (^|) the commutator 
L+L - LL + = L\(t)^og\u(x,t)\ 2 ] xx = L\{t)(h - h x ) = -L\{t)A 01 (x,t). So we have 
g^-g^=iL\{t)A^t). 

If the transformation function u is nodeless, the potential difference (|5|) has no 
severe singularities and the operators L (^) and L + (|?p are well defined. There exists a 
sole (up to a constant factor) nodeless function in the space H°(R). Beyond this space 
there are many nodeless functions suitable for use as transformation function. 

Let us assume that the nodeless transformation function u is such that the absolute 
value of m _1 (x, t) is square integrable. Then the solution of Eq. (§) of the form 
(HP is the discrete spectrum eigenfunction of In the space 7\ the operator L 

induces a subspace H n (R) = {ip : tp = Lip, Vif> G H°(R)}. If we denote H W (R) = 
kerL + = {jv, V7 G C} then we have the decomposition of the Hilbert space H 1 (R) = 
H 10 (R) © H U (R) C Tj. In this case the lowest eigenvalue of the super-operator 5* 
is nondegenerate and its eigenfunction ^/q = N (0, v) T (N being the normalization 
constant) represents the vacuum state and is annihilated by both supercharges Q, Q + 
and, consequently, we have the exact supersymmetry. 

We can repeat all these reasonings starting from the Eq. (0), Hamiltonian hi and 
symmetry operator g^ 1 ' and construct new Hamiltonian fi2 = hi + Ai^(x, t), symmetry 
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operator g^> = g^- 1 ' — iL\ 2 (t)Ai t 2{x,t), and new exactly solvable Schrodinger equation 
etc. The chain of N time-dependent Darboux transformations can serve us as a basis 
for construction of iV-extended supersymmetric model or parasupersymmetric one. 
These constructions can be made in full analogy with the corresponding stationary 
constructions 0-0, fll3 . 

4. Let we have a chain of the well defined Hamiltonians ho — > hi — > . . . — > ob- 
tained by the subsequent applications Darboux transformation operators L Q i, L 12 , . . . , 
Ln-i,n- With the help of Eq. (|3]) we can express all the intermediate transformation 
functions through the solutions of the initial equation (JT|) . As a result, we eliminate 
all the intermediate Hamiltonians and Schrodinger equations and pass from Eq. (JJ) 
immediately to the final one. The chain of the first-order transformations collapses in 
a single iV-order transformation and we obtain a nonstationary analogue of the well 
known Crum-Krein formula 



14 , 15 



Lo. 



N 



L N -i,nL N - 2 ,n-i ■ ■ ■ L 0t i = L N (t)W(ui, u 2 , ■ ■ ■ , u N ) 



Ui 


U 2 ■ 
U 2x ■ 


1 

• d x 


u lx 


u 2x 


. di N) 



where the conventional symbol W(u\, u 2 , . . . ,un) stands for the Wronskian of the trans- 
formation functions Mi, m 2 , . . . , un- If there is no necessity to introduce the intermediate 
Hamiltonians and Schrodinger equations we can avoid the reality condition for every 
transformation function Ui and impose a single condition of reality of the final potential 



log 



/ W(u 1 ,u 2 , ...,%) 



\W*(u 1 ,u 2 , ■ ■ .,u N ) 
The potential difference for this case reads as follows 

A ,n = ~ [log \W(ui,u 2 , ■ ■ -,u N )\ 2 ] 
and function L^{t) can be chosen real 



0. 



11^ 



L N (t) = exp I 2 dtlm [log W(ui, u 2 , ■ ■ ■ , u N )} 



We want to stress that the condition (y) may be violated but the condition ( pTT|) should 
hold. The analogous case for the stationary Schrodinger equation is called the irre- 
ducible one 0] and the operator Lo.v is said to be the Nth-order Darboux transformation 
operator JT^ . 



5. Consider the case when all the transformation functions 
eigenfunctions of the same symmetry operator g(°\ ig^Uk 



1,2, 



iV are 



P 



1,2, 



and every Wronskian W(ui,u 2 , 
in R for all t and satisfies the reality condition 



a k u k , ct x > a 2 > ■ ■ ■ > a N 
N, W{ui) = U\ conserves its sign 
). Together with the chain of the 



5 



Schrodinger equations with the well-defined Hamiltonians hi we now have the chain of 
the well-defined symmetry operators g^K We call this case the completely reducible one. 
Each pair of the adjacent operators and g( p+ V as well as the Schrodinger equations 
with the adjacent Hamiltonians h p and h p+ i is intertwined by the well-defined first- 
order transformation operator L P)P +i, p = 0, 1, . . . , N — 1. Consider the nth-order 
(n = 1, 2, . . . , N) operators L p , p+ri = L p+ri _i :P+n . . . L p+liP+2 L PjP+1 . It is evident that 
every operator L PtP+n intertwines two Schrodinger equations with the Hamiltonians 
hp, h p+n and two symmetry operators g^ p \ g( p+n ) and consequently is the nth-order 
Darboux transformation operator. If we denote the operator Laplace adjoint to L PtP+n 
by L p4>+n then we can establish the following factorization properties: 

n 

5=1 • +n (12) 
Lp : p-\- n L ptP+n = Y[ (W ^ ^ ~ CYp+k) ■ 

k=l 

The whole chain of the transformations is associated with supercharges Q p , q = 
L p , q e Ptq , p < q; p,q = 0,1,..., N, super-Hamiltonian H = diag{/io, hi, . . . , h??} and 
super-operator S = dia,g{ig^°\ ig^\ . . . , ig^}. The notations e M (see ||) refer to 
(N + l)-dimensional square matrix whose rows and columns are labeled from to N 
and contain zeros everywhere except units at the intersection of the pth column and 
qth row. As Q^ q we denote the adjoint operators Qp q = L^ q e q>p . The chain of the 



iV + 1 Schrodinger equations is summarized in a single matrix equation (|10f) where 
/ is the (N + l)-dimensional unit matrix. The intertwining relations between the 
operators L PjP+n , L£ +n and the Schrodinger operators id t — h p , id t — h p+n translate into 
the commutation relations between Q PtP+n and matrix Schrodinger operator ild t — H. 
Hence, the supercharges Q PtQ and are integrals of motion for the system with the 
super-Hamiltonian H and the symmetry operator S commutes with all Q p , q s and Q pq s. 
The reducibility condition leads to the nonlinear algebra of the following relations: 

Q P , q Qs,p = Q s ,q, s<p<q, s,p,q = 0,l,...N, 

n 

Qp,p+nQp t p+n+m Yl. W a p+i)Q p+n,p+n+mi P n m — 
1=1 

n 

Qp—n—m,pQp—n,p ^p+i— l)Qp— n— m,p— ni P n m — 0; P — -^"> (13) 

i=l 

n 

Qp,p+nQp.p+ n Qp,p+n = Yl ~~ a p+i)Qp,p+n, p + U < N, 

i=l 

n, m — 1, 2, . . . 

and of the Hermitian conjugated ones. All the other products of every two supercharges 
are equal to zero. Hence we obtain the nonlinear algebra with a parasupersymmetric 
structure. 

6. As an example consider the free particle parasupersymmetric model: Vq(x, t) = 
0, R — R. Functions ip n (x,t) of the discrete basis set in the space H°(R) are eigen- 
functions of the operator |T2[: gW = K_ 2 - K 2 = i{l + t 2 )d 2 x + txd x + t/2 - ix 2 /A, 
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g^ipx = i^ipx and correspond to A = A„ = — n — 1/2, n = 0, 1, 2, . . . (■?/>„ = Vx)- Their 



coordinate representation is as follows (12 



V>aOM) = (1 + t 2 )" 1/4 exp[zx 2 t/(4 + 4t 2 ) +iAarctant]g A (^), 

^ = x/ v / m 2 " ^ ^ 

where Q\{z) is the parabolic cylinder function satisfying the equation 

Ql(z)-(z 2 /A + \)Q x (z)=0, 

where prime denotes the derivative with respect to z. It is easy to see that the functions 
fllil) for any real A satisfy the reality condition (D). 



Functions ([14]) at A = A m = m + 1/2, m — 0, 1, 2, . . . do not belong to the space 
H°(R) and at even m are nodeless for all i e M and suitable for use for the single 
Darboux transformation as transformation functions u m = ip\ m - (Note that for odd 
m they are equal to zero at x — 0). Using formulas (043) we find the transformation 
operator 



L = L ,i = Li(t)<9 x - - W- - «m— — t-^, = VI + t 2 

2 V 1 — it He m {iz) 

where ife m (z) = 2" m / 2 J ff m (z/ v / 2) (# m (z) are Hermite polynomials). New Schrodinger 
equation potential is calculated by the Eq. (|5|) 

m—2{i%) n 2~^^m— 1(^) -n 

^ = ^ = (1 + [M--l)^^-2m^- r -l]. 

The same potential difference defines the symmetry operator g( l > = g(°> — iLl(t)A 0jl 
of the Schrodinger equation with the potential V\. Its discrete eigenfunctions obtained 
with the help of the operator L ,i and with the use of formula (|) form the basis of the 
space H l (R) of the square integrable solutions of the new Schrodinger equation 



<p = v / ^!(27r)- 1 /4 [jLl(t)<(X)t)] -i ) 

# (1 Vo = i(m+ 1/2)^0, 
¥Wi = [™K n + ™ + + it)v / 27r]~ 1/2 exp[-x 2 /(4 + 4it) - marctant] 
x[He n+1 (z)+imHe n (z)^^l 
# (1 Vn+i = -i{n + l/2)<f n+h n = 0, 1, 2, . . . . 

The normalisation constants of the functions LQ^if) n , n = 0, 1, . . . are found with the 
help of the following relatrions: 

(L ,iVr» I AuVU = (VVi I ^(tl L 0,l | V'n) 
= (^n I ^ (0) + m + 1/2 | ^n) = (n + m + l)(ip n I ^n)- 

The normalisation constant for the function tpo is obtained by the direct calculation. 
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It can be shown that the action of the operator Lq,i on the functions ([14]) with 
A = Xi = I + 1/2, I = m + 1, m + 3, . . . gives nodeless solutions of the new Schrodinger 
equation 

vi(x,t) = L 0>1 ip Xl (x,t) = (1 -zt)- 1 / 2 exp[x 2 /(4-4it) + z/arctant]^gy, 
fm,l( z ) = i[Hei(iz)He m+ i(iz) - He m (iz)Hei +1 (iz)} 

which do not belong to the space H 1 (R) and are suitable for use for the second trans- 
formation as the transformation functions. Note that the absolute value of the function 
(x, t) is square integrable for all t and I = m + 1, m + 3, .... This fact can easily be 
shown by the direct calculation of their normalisation constants. The transformation 
operator Li i2 constructed with the help of the transformation function u(x, t) = Vi(x, t) 
according with the formulas ([3] - f|) has the form 



L\ 2 = VTTT 2 d x - ~ A /^ + im^4$ - i^ {z) 



2 V 1 - it He m (iz) f m i(z)' 

New potential difference A = A 1>2 = V 2 — V\ is anew calculated by the Eq. (|5|) and 
leads to a new exactly solvable potential 



T/ o/l I j2\-l\-, I fml( Z ) ( f'ml( Z ) \ 



New Schrodinger equation symmetry operator has the form = — iL\{t)Ai^ = 
g (0) - %L\ (t)A 0t2 = 5 ,(0) - iL\(t)V 2 - (Note that since ui/u* = Vi/v*, the time-dependent 
factor Li(t) calculated by the formula (|4]) is always the same for this case). Its or- 
thonormal eigenf unctions Xi+i( x ->t), i = 0, 1,2, . . . forming together with the function 
Xo{x,t) ~ [Li(t)v*(x, t)]^ 1 the discrete basis set of the space H 2 (R) of the square inte- 
grable solutions of the Schrodinger equation with the potential V% should be found by 
applying the operator Li j2 on the functions ifi(x,t). Finally we obtain 



2 



X0=[(27T)- 1 /2/ ! (/_ m )(l_, t )-l]l/^ 

— x J (4 + Ait) + il arctan t\ - 



fml(z) 

Xl = [(2vr)- 1 / 2 m!(/ - m)(l - it)' 1 ] 1 / 2 

r 9 u a a ■ \ ■ ,Hei(iz 
x exp| — x J [A + Ait) + im arctan tj— — 



fml(z) 

Xn+2 = [(n + I + l)(n + m + 1)]~ 1/2 
x [-(I + n + l) Vn + m)u m Z { ^ Ul \ ] 
ra = 0,l,2,..., 



where 



W{(p n ,ui) = (n\V2n)- 1/2 (l + t 2 )' 1 exp[itx 2 / (2 + 2t 2 ) n) arctan t] x 

[nHe n -x(z)Hei{iz) + iHe n (z)Hei +1 (iz)] 
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and 

W(u m ,ui)= ((1 - it)Vl + t 2 ) f ml (z) exp[x 2 / (2 - 2it) + i(m + Z) arctant]. 

Second-order operator L ,2 = ^1,2^0,1 intertwines the initial Schrodinger equation 
(with the zero potential) and the final one (with the potential V2) and participates in 
the following factorizations [see Eq. ©]: Lj 2 L ,2 = (i# (0) + m + 1/2) (i^^ + I + 1/2) 
and L 0) 2-^o 2 — (w® + m + 

l/2)(ic/( 2 ) + / + l/2). The operators L 0: 1, L12 and L 0) 2 form 
the complete set of the transformation operators for the case N = 2. We associate 
with this set the complete set of the supercharges Q PtQ = L Pjq e Pyq and Q+ q = (Q p , g )^ , 
p, q — 0, 1, 2 (p < q), which together with the super-operator 5 = i di&g{g(°\ g^} 
form the following nonlinear algebra [see formulae fll3|) 1: 

[£> Qp.g] = 0) <5i,2<5o,i = Qo,2, 

Qo.iQoiQo.i = (5 , + m + l/2)Q 0) i ) 

Ql,2Q^ 2 Ql,2 = (S + / + l/2)Ql,2, 

QoM 2 Qo,2 = (S + m + 1/2)(S + I + 1/2)Q ,2, 
<3o,iQo;2 = (5 + m + l/2)Q+ 2 , 
Q{ 2 Qi )2 = (S + Z + l/2)Q+ 1; 

and relations Hermitian conjugated to these. All the other products of any two super- 
charges are equal to zero. 

The lowest eigenvalue —I — 1/2 of the superoperator S is nondegenerate. Its eigen- 
function (vacuum state) has the form ^ = (0, 0, Xo) T ■ Two lowest excited states of this 
operator = (0, 0, Xi) T , ^1 = (0, <po, 0) T ] correspond to the eigenvalue — m — 1/2. 
All the other eigenvalues are triple degenerate and equal n + 1/2, n = 0,1,2,.... 
Their eigenf unctions have the form ^> l n+2 = (0, 0, Xn+2) T , ^l +2 = (0, <f n +i, 0) T , \E^ +2 = 
(V>«,0,0) T . 
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